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Abstract 

The mean electromotive force (MEMF) in a rotating stratified magnetohydrodynamical 
turbulence is studied. Our study rests on the mean-field magnetohydrodynamics framework 
and T approximation. We compute the effects of the large-scale magnetic fields (LSMF), global 
rotation and large-scale shear flow on the different parts of the MEMF (such as a - effect, 
turbulent diffusion, turbulent transport, etc.) in an explicit form. The influence of the helical 
magnetic fluctuations which stem from the small-scale dynamo is taken into account, as well. 
In the paper, we derive the equation governing the current helicity evolution. It is shown that 
the joint effect of the differential rotation and magnetic fluctuations in the stratifled media 
can be responsible for the generation, maintenance and redistribution of the current helicity. 
The implication of the obtained results to astrophysical dynamos is considered as well. 

1 Introduction 

The mean- field magnetohydrodynamics presents one of the most powerfu l tools for explorin; 
the nature of the l arge-s cale magnetic activity in cosmic bodies ( Moffatt . 1978t Parkeil 197' 



Krause and Radler . 1980l ). It is widely believed that magnetic field generation there is governed 



by interplay between turbulent motions of electrically conductive fluids and global rotation. The 
growth and evolution of the large-scale magnetic fields (LSMF) in cosmic plasma strongly depends 
on the mean electromotive force, £ — (u x b) , which is given by the correlation between the 
fluctuating components of the velocity field of plasma, u, and the fluctuating magnetic fields, b. 

The global rotation, stratification and the strong LSMF can substantially modify the structure 
and amplitude of the mean electromotive force (hereafter, MEMF) leading to the rich variety 
of the turbulence effects driving the evolution of the LSMF in cosmic bodies, e .g., the a-effect 



( Roberts and Sowardl . [l975l : lMoffattl . ll978l:lKrause and Radleil . ll980HParkeiill979l : lRudiger and Kichatinov 

the rotational l y-induced anisotropy of turbulent diffusion and e ffectiv e drift of LSMF 
([Roberts and Soward . 1975t Krause and Radler . 1980l : Kichatinov et al. . 1994), etc.. Generally 
speaking, the nonlinear effects of the small-scale Lorentz forces on the MEMF and LSMF evolu- 
tion stem from two sources. One is driven by perturbations of the LSMF due to turbulent motions 
and another one is due to magnetic fluctuations, which are maintained by the small-scale dynamo 
action in a turbulent medium. The role of the small-scale dynamo in the LSMF evolution is still in- 
sufficiently understood^Ni^^ 

ture, e.g., (|Moffattl .l Il978t iFrisch et al.l . Il975l : IPouauet et al.l . Il975l : [Brandenburg and Subramanianl . 
20051 ) . According to the mentioned studies the most important effect of the growing magnetic fluc- 
tuations on the LSMF evolution is caused by the helical part of magnetic fluctuations. The mag- 
netic helicity conservation law, if applied to the mean-field magnetohydrodynamics, requires that 
the amount of helicity contained in the LSMF (controlled mostly by a-effect) should be roughly 
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the s ame and opposite in sign to its counterpart in the smaU scales, see (jKleeorin and Ruzmaikid . 
19821 ). In this way the helical part of magnetic fluctuations, which is excited both due to shred- 
ding the LSMF by turbulent motions and due to small-scale dyri a mo, e ff ectively saturates th e 
generation of the LSMF by a-effect (Vainshtein and Kitchatinovl (1983); Brandenburg ( 2001 ): 
iPield and BlackmanI ( 20021 ) : lBlackman and Brandenburg! ( 20021 )'). Further discussions on this sub- 
ject can be found in above cited papers. Their main lesson is that the construction of the realistic 
mean-field dynamo theory requires the evolution of the small scale magnetic (or current-) helicity 
to be taken into account. 

Currently, there are two basic schemes for computing the MEMF of turbulent fields. One is the 
quasi-linear approximation (the same approximation is called the FOSA or SOCA in literature). 
A com prehensive discussion about its applicability and val i dity in astrophysics can be found at pa- 
pers bvl MoffafrO ( 1978[ ): Parker! ( 19791 ): Krause and Radler ( 1980l ): Brandenburg and Subramanian! 
(!2005). This scheme remains one of the main tools of the mean- field magnetohydrodynamics. How- 
ever, one of unfortunate problem of SOCA is that the contribution of the magnetic fluctuations 
(and the corresponding magnetic helicity) driven by the small-scale dynamo is hardly possible to in- 

clude in the th e ory in self-consistent way. The third order closure s chein e based on r-approximation 

Orszag! ( 1970! ): Vainshtein and Kitchatinov! ( 1983! ): lRadler et aD ( 20031 ) : lBrandenburg and Subramanian 



mm) 

gi ves a chance to consider, roughly, the effects of the small-scale dynamo on the MEMF. 



Following ( Brandenburg and Subramanian . 2005! ) (hereafter BS05), I will call it MTA ( minimal tau 
approximation). Different kinds of this approximation are used in the literature, see dVainshtein , 



1983 [ Vainshtein and Kit chatinov , 1983 [ Radler et al.!. 120031: iBrandenburg and Subramanian! . 12005 ; 
Rogachevskii and Kleeor in. 2003: iBlackman and Field . 2002! : Field and Blackman . 2002). In the 
paper we follow procedure described in BS05. Furthermore, the variant of tau approximation 
with a scale-independent relaxation time, r, is applied. For thi s reason, some results obtai i ied in 
the pa p er can be different o f those that are given elsewhere: (jRogachevskii and Kleeorinl 12003 
2004bl la! 



Radler et all . 120031 ) 



The main purpose of this paper is to compute the MEMF via MTA taking into account the 
influence of the global rotation and LSMF on the turbulence. The stratification of the medium 
and the large-scale shear are taken into account as well. The influence of rotation, LSMF and 
uniform shear on the different parts of the MEMF (such as a - effect, turbulent diffusion, turbulent 
transport and etc.) is explicitly defined via factors describing the efficiency of rotational and LSMF 
feedback on the turbulent flows. The influence of rotation is measured by the Coriolis number, 
f2* = 2Qtc , where Q is the solid body angular velocity and Tc - the typical correlation time of 
turbulent flows. The influence of LSMF is measured hy (3 — B/ (^Ucy/]Tp), where B is the strength 
of the LSMF, Uc is a typical rms velocity of turbulent flows and /i, p are the magnetic permeability 
and the density of the media, respectively. Following the basic approach developed in above cited 
papers we derive the equations governing the evolution of the current helicity both in rotating and 
in magnetized turbulent flows with imposed uniform shear. 

The paper is structured as follows. In the next section we shortly outline the basic equations, 
assumptions and the computational scheme for derivation of the MEMF and the evolutionary 
equation for current helicity. Section 3 is devoted to the results of calculations of the MEMF 
for different situations (slow rotation, strong LSMF, vice versa and etc.). In section 4 we derive 
the evolutionary equation for current helicity. In section 5 we summarize the main results of the 
paper. 



2 Basic equations 

In the spirit of the mean- field magnetohydrodynamics, we split the physical quantities of the 
turbulent conducting fluid into the mean and randomly fluctuating part with the mean part 
defined as the ensemble average of the random fields. One assumes the validity of the Reynolds 
rules. The magnetic field B and velocity of motions V are decomposed as follows: B = B + b, 
V = V-t-u. Hereafter, everywhere, we use the small letters for the fluctuating part of the fields and 
capital letters with a bar above for the mean fields. The angle brackets are used for the ensemble 
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average of products. Foll owing the lines of two-scale approximation ( Roberts and SowardL Il975t 
Krause and RadleJ . llQSOl ) we assume that the mean fields vary over the much larger scales (both 



in time and in space) than the fluctuating fields. The average effect of the MHD-turbulence on the 
LSMF evolution is described by the MEMF, £ = (u x b). The governing equations for fluctuating 
magnetic field and velocity are written in a rotating coordinate system as follows 



drrii 



dh 

'dt 

2 (J7 X m),. 



= V X (ux B + Vx b) +?7V2b + (S, (1) 

j/Amj + 1/ (G • V) TOi (2) 



/ 2 , (b-B) 



+ -Vj {Bjbi + B,bj) - Vj {Vjm, + V^nij) + /, + g-,, 



where 0,5^ are nonlinear contributions of fluctuating fields, m = pu, G = Vlogp is the density 
stratification scale of the media, p - the fluctuating pressure, O, - the angular velocity responsible 
for the Coriolis force, V - mean flow which is a weakly variable in space, f - the random force 
driving the turbulence. 

To compute £ it is convenient to write equations U]) and ^ in Fourier space: 



d 



B, 



B, 



%(z - q) j (q) - miiz - q) \^-^ j (q) 
6i(z - q)t,(q) - 6,(z - q)t,(q)l dq + 



dq (3) 



d 

— + i^z^ + iiy (Gz) ) TO, = (Oz) (z X m) . 



(4) 



mif{z)zi / mi(z - q)V"/(q) + TO/(z - q)yi(q) 



dq 



+ / bi{z-q)Bf{q) + bf{z-q)Bi{q) dq, 

where the turbulent pressure was excluded from by convolution with the projection tensor 
7r.y (z) = Sij — ZiZj, Sij is the Kronecker symbol and z is a unit wave vector. The equations for 
the second-order moments which make contributions to the MEMF can be found from (|3l4p . As 
the preliminary step we write the equations for the second-order products of the fluctuating fields, 
and make the ensemble averaging of them, 



d I 

— im, (z) hj (z') 



r/i,^ (z,z') - (?7z'2 + vz^ + \v (Gz)) (m, (z) (z') 



(5) 



{rrii (z) TOj(z' - q)) 
(m,(z)TO,(z'-q)) (^)(q) 



(q)- 



dq - 2 (Jlz) e^inZi (?7i„(z)6j(z') 



TO, (z) h (z' - q) ) V, (q) - ( rn, (z) 5, (z' - q) ) 1/, (q; 



dq 



^iliz - q)fe, (z') ) \^/(q) + (%(z - q)6, (z')) ^^Kq) dq 

dq, 



- \'Kij{z)zi 

+ -^zi^,s(z) I {hi{z - q)fo, (z') ) Bf (q) + ( &/(z - q)6, (z') ) Bi {q) 
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dt 



(m^ (z) TOj (z')) = -2 (r2z) £i;„f; (to„(z)toj(z')) - 2 (fiz') ej;„Z;' (mj(z)?7i„(z')) 

- i7i"»/(z)2; / (m,(z-q)mj(z')> V/(q) + (m/(z-q)TOj(z'))V^i(q) 



(6) 



dq 



i7r,-f z z, 



(z) m,(z - q)) Vf{(\) + (mj (z) r%(z - q)) y,(q) dq 



+ -^■^f{^)zil {bi{z~(i)mj{z'))Bf{q) + {bf{z-q)mj(z'))Bi{q) 



dq 



mi (z) 6; (z - q) ) 5/ (q) + ( rhi (z) 6/ (z - q) ) (q) dq 



+ T/i^,(z, z') - 1/ + + z (Gz) + t (Gz')) (to, (z) m, (z')) 
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S.(z)6,(z')) = Tft^.(z,z')-(ryz'2+r7z2)(&.(z)6,(z')) 



5/ 



6, (z) m,(z' - q) ^ (q) - 6, (z) m^z' - q) ^ (q) 



5, 



m,(z - q)6, (z') ) ( ) (q) - (™,(z - q)6,. (z') ) ( ^ ) (q) 



(7) 
dq 

dq 



6, (z) 6Kz' - q)) V,{q) - (^h (z) b,{z' - q)) V^q) 



dq 



bi (z - q)S,(z')) ^^.(q) -{b,{z- q) 6,(z') ) ^/(q) dq, 



where, the terms xh\^'^'^^ involve the third-order moments of fluctuating fields and second-order 
moments of them with the forcing term. 

To proceed further, it is convenient to introduce some notations which are used in the literature. 
The double Fourier transformation of an ensemble average of two fluctuating quantities, say / and 
(7, taken at equal times and at the different positions x, x', is given by 



(/(x)5(x')> 



/ (z) g [z')) e'(^-^+^'-''')d3zd3z' 



(8) 



Let us define the "fast" spatial variable r by the relative difference of x,x' coordinates, r = x — x'. 
The "slow" spatial variable R is R = (x + x') /2. Then, eq. ([5]) can be written in the form 



(/(x)5(x')) = 



1_ 



1 



/ k + -K .g -k + -K ) c'(^-^+'' ""M3Kd3k 



(9) 



where I have introduced two wave vectors: k = (z — z') /2 and K = z -I- z'. Following BS05, we 
define the correlation function of f and g obtained from ([9|) by integration with respect to K, 



$ (/, g, k, r) = I l^f (^k + Ik) g (^-k + Ik) ^ e'(^-^)d^K. 



(10) 



For further convenience we define the second order correlations of momentum density, magnetic 
fluctuations and the cross-correlations of momentum and magnetic fluctuations via 



f)y(k,R) = $(TO„TOj,k,R),p2(u2^ (R) = y"{),^ (k,R)d3k, 

%(k,R) = $(6„6j,k,R),(62)(R) = y" /i,, (k, R) d^k, 

kij (k, R) = $ {rhi , 6^ , k, R) , pS^ (R) = (k, R) d^k. 



(11) 
(12) 

(13) 
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Let us now return to equations 



and (l7|). As the first step, we approximate the Th. 



terms by the corresponding r relaxation terms of the second-order contributions, 



Th 



(>') 



Th. 



(v) 



Th 



(h) 



rhi (z) hj (z')^ /tc, 
rrii (z) rhj (z')) - {rhj (z) rh.j (z'))^°^ 



k (z)&,(z'))-(6,;(z)&,(z') 



(0) 



(14) 
(15) 

(16) 



where the superscript denotes the moments of the background turbulence. Here, Tc is indepen- 
dent on k and it is independent on the mean fields as well. Furthermore, for the sake of simplicity, 
we restrict ourselves to the high Reynolds numbers limit and discard the microscopic diffusion 
terms. As the next step we make the Taylor expansion with respect to the "slow" variables and 
take the Fourier transformation, pl}|) . about them. The details of this procedure can be found 
in BS05. In result, we obtain equations for the second order correlations of momentum density, 
magnetic fluctuations and the cross-correlations of momentum and magnetic fluctuations. 



dk, 



dt 



-i (Bk) y-f-^ 

-•■JlVil-Dj H 

P M 



2 



BV) Iv^^h 
kiBij d 



{Bkl^dv,, (GB) 



2p 



2p 



dkf 
dk. 



kikfBfihij 



2^ (nkj kpSipi (k^^ >fij + ^£ipi 



dkf Tc 
ilfc j S/pkij + kp (ilV) kij j , 



(17) 



dt 



.(0) 

-2 Ink) kp {Sipivij + ejpiVii) ^ "^n^jJ 



+ 2kfVf^i [kiVi] + kjVii^ + kiVfj 



dkf 



i(Bk) iki. 



1 



Bi {kij,i -f k*i i) + B^^ik*i + Bj^ik^^i ~ 2kfBf^i {kik*i + kjkij 

kp {{vlv) - 2 [vtk^ (i^^)) 



2 

Bij 



2 



+ (ilfc Vp {eipivij - EjpiVii) , 



(18) 



% = + + (19) 

\{BV) (BG)\ (BA fBA . 



2p 2p I ^ V P 

l'Bi\ d{k,j 
ki- 



2 \ p J J.- dkf 

where k*^ — $(6j, m^, k, R), k is the unit wave vector, the indexes behind the comma stand for 
th e spatial derivatives. Equati ons (I17I18I19P are in agreement with those considered in the paper 
by iRogachevskii and KleeorinI ([2004a,). 
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To solve (|17|18|19p we neglect the time derivatives at the left hand side of equations and apply 
the perturbation method. The mean field inhomogeneities and stratification scales of turbulence 
are considered as small. We shall not reproduc e exphcitly th e rather bulky derivations which are 
explained elsewhere: Rogachevskii & Kleeorin(l2003l l2004bl) . The solution of (|17ll8ll9|) will be 
given for two specific cases. In the first case we apply no restriction to the angular velocity ( 
the Coriolis number, 17* = 2Utc, is arbitrary) and LSMF is assumed to be weak. In the second 
case we keep the linear terms in angular velocity and solve eqs. (|I7|I8|f 9p for the case of arbitrary 
P = B/ (uc^/Jip), where B is the strength of the LSMF. In all derivations we keep contributions 
which are the first order in the shear. Furthermore, for the contributions involving the shear we 
make two additional simplifications. The first one is that we neglect the density stratification, 
but leave the contributions of the turbulence intensity stratification. Additionally, we discard the 
joint effect of the Coriolis force and the shear to the MEMF. In the present study I consider an 
intermediate nonlinearity which implies that effect of the mean magnetic field and global rotation 
is not enough strong in order to affect the correlation time of turbulent vel ocity field. 



F or integration in k-space I adop t the quasi-isotropic form of the spectra (jRoberts and Soward 



19751 : iRtidiger and Kichatinovlll993( ) for the background turbulence. Additionally, the background 
magnetic fluctuations are helical, while there is no prescribed kinetic helicity in the background 
turbulence: 



-,(0) 



;(o) _ 



Ti-y (k) + ^ ihVj - kjVi) 



8^ ' 
g(fc,R) 



. AA(fc,R) 



(20) 
(21) 



where, the spectral functions E{k,H), B{k,'El),J\f{k, R) define, respectively, the intensity of the 
velocity fluctuations, the intensity of the magnetic fluctuations and amount of current helicity in 
the background turbulence. They are defined via 



where h. 



,(0)2 



(0) 



Ank^ 



g(fc,R) 



d-'k, K 



(0) 

c 



1 

HP 



:^d3k, 

47rfc2 



(22) 



— ■ ■ V ^ u ■ I / {pp) . In final results we use the relation between intensities of 

magnetic and kinetic fluctuations which is defined via i3(fc,R) = eppE {k,'R). The state with 
e = 1 means equipartition between energies of magnetic an kinetic fluctuations in the background 
turbulence. The point to note is that inconsistency between ([^0]) and (^1]) does not influence the 
final results. The general structure of the mean electromotive force vector obtained within the 



given framework are in agr e emen t with the known results from the literature (jRadler et al.l . 12003 



Rogachevskii and KleeorinL 2003 ). We keep the current helicity contribution in the background 



turbulence to investigate the nonlinear saturation phase of the helical large-sca le dynamo. 

Th e final remarks in this section concern with discussion given in the paper bv lRadler and RheinhardO 
(|2006l ). There, authors argue that r approximation may lead to results which are in conflict with 
those of SOCA. One difference is apparent between the two approaches: there is no overlap in ap- 
plicability limits of SOCA and r approximation. The given scheme to obtain ()17|18|19|) is hardly 
justified for small hydrodynamic Reynolds numbers. The same is true in a highly conductivity 
limit, where SOCA can be valid only for the small Strouhal numbers. Currently, the range of r 
approximation validity is purely understood. This problem requires further careful study. 

There is another reason for difference between results presented in the paper and those of 
SOCA. In the given variant of r approximation the relaxation time Tc is independent of k. This 
issue is especially important in computing effects of the nonuniform LSMF and shear. Perhaps , 



the spectral r-approximation can correct this defect. For more detail, see (jRadler et al.l . 12003 



[Rogachevskii and Kleeorin . 2007 ). Hense, in confronting MTA and SOCA, it is of some use to 
simplify the expressions obtained within SOCA by applying the mixing-length approximation. 
The transition from SOCA to MLT can be done by replacing the spectrum of turbulent fields by 



the single-scaled function of the form 6 {k 



S (uj), and applying rjk = vk = , here ic 



IS 
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the c orrelation length of the turbulence. For more details, see (jKichatinovLll991tlKichatinov et all 
1994h . 



3 Results 



3.1 Weak LSMF, arbitrary Coriolis number 
3.1.1 Spatially uniform LSMF 

We divide the electromotive force into different contributions, in particular, S^""^ contains the 
effects of stratification, and is due to shear. The contributions due to shear are computed 
only in slow rotation limit. We find the following expression for S'-"-^: 



= {(£-l)(/('^)(UxB)+/('^)(e.B)(exU))+/i")(GxB) 

+ /l"^) ((e • G) (e X B) + (£ - 2) (e • B) (e x G)) 

+ ft^e (e . B) (e . U) + /^^'b (e • U) + f^^e (e • B) (e • G) 

+ ft^ (e (B . U) + U (e . B)) + ft^ (e (B • G) + G (e • B)) + ^"^B (e • G) 



(23) 



+ /('^'(e(B.U)-U(e.B))+/('^)(e(B.G)-G(e.B))}( 
+ 2{/('')B-/('^)e(e.B)}r,4"\ 

where functions /|°| = /|°j(ri*,e) (and all which are used below) are given in Appendix A, 

U = V log (it*^")^) is a scale of the turbulence intensity stratification, e = ll/|f2| is a unit vector 
in direction of global rotation. For the slow rotation limit {il* 0) we get : 



f^-^^ln-^o = aoB+(u 



,(0)2 



(e-1) 



U X B) 



(G X B 



(24) 



,(0)2 



r,— {(£ + 2) (((G X e) X B)) + (e + 1) (((U x e) x B)) } 



2e((e-U) + (e-G))r2* 2{e-V)Q* 4 (e • G) 17* 




15 



l(0) 

3 



(e,Gj + e^Gi) (e + 4) + {e,Uj + CjUi) e + — 
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(25) 



are 



where only linear terms in $7 are kept. Except c ontri butions due to G equations (1241) and \ 
in agreement with results bv lRadler et al. ( 2003 ) and Brandenburg and Subramanian ( 2005[ ). The 
mean transport of the LSMF due to stratification of turbulence is give n by second term i n ((2^ . 
They are in agreement with the mixing-length expressions obtained bv iKichatinov ( 1991 ). Note 
that, additional components of the turbulent transport may be excited due to the antisymmetric 
part of a-tensor in ([25]) . 

For the fast rotation limit {Vl* oo) of (P5)l we get 



5(a) I 




,(0)2 



(e-U) 



+ (e-G 




B)) 



(26) 



where, we keep the next order contribution in U,* for the current helicity, as well. The reason for 
this will be clarified later in section [H Except the helicity term, eq.(l26l) is in identical agr eement 
with the mixing- length approximation results obtained bv lRiidiger and Kichatinov (1993) within 
SOCA. 

In the case of the spatially uniform LSMF the shear contributions to the mean electromotive 
force are expressed as follows: 
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Figure 1: The modification of standard alpha effect (cf. iKrause and Radlei ( 1980? )) due to shear. 
The helical motions (denoted with (u • V x u)) go up, drag and twist the LSMF B^, where index t 
denotes the toroidal component of LSMF. The shear, AVp, additionally, folds the loop in direction 
of large-scale flow. The effect is equivalent to inducing the transversal large-scale electromotive 
force, £p (here index p denotes the poloidal component of the MEMF), and the magnetic field, 
parallel to original one. Direction of the induced field depends on the sign of the helicity. For 
the situation given on the picture, the induced field B^ quenches the original LSMF in direction 
of the gradient of the mean flow. This means that the LSMF is effectively pumped in opposite 
direction. 



fi'^ = Emm {A4UkBnV^^k+A2BkVn,kUm+A3 (B • U) V rn,n + A^V k,nB kU,n} (27) 

! (0) -, o 

+ r!^^ (W X B)^ - ^^r!hf {F„, +F,„} B„, 

where W = V x V, we assume that (U • V)V = and Ai = (2e - 1) t?/15 Ai = - (Se 1)t.V15, 



A3 = {e -I- 1) /6, = —A3,. Coefficients Ai_3 correspond to those from lRiidiger and Kichatinov 



( 20061 ) (hereafter RK06) and A4 is corresponding to th eir A5. Recently, simil a r con tributions of 
the large-scale shear were calculated within SOCA by iRadler and Stepanov ( 20061 ) (RS06), as 
well. We have to note that both the RK06 and RS06 results are related with the case £ = 0. 
The ([27| differs with results obtained in RK06 and RS06 papers. For example, after applying 
the mixing-length relations r\y? ~ vk^ = to expressions given by RK06 we get Ai = rf/3 
(in our case — t^/15) and A2 ~ — t^/60 (compare to our — t^/15). Unfortunately RK06 did not 
present the results for other coefficients. The comparison with RS06 is given in Appendix B. The 
difference between the given results and those by RK06 and RS06 can be explained, in part, by the 
crudeness of the given version of tau approximation. Here, we assume that Tc is independent of k. 
This especially influences the accuracy of calculations of the contributions due to shear because 
they involve the derivatives in k space. 

According to ((77|) the joint effect of current helicity and shear contributes to pumping of LSMF. 
The interpretation of the effect is difficult to illustrate. To show the general idea we invoke an 
auxiliary illustration of effect for the helical turbulent motions. It is shown on Fig[TJ 



3.1.2 Anisotropic diffusion, the O x J and shear-current effects 



In rotating turbulence the magnetic diffusivity become anisotropic (jKichatinov et al.l . 119941 ). The 
corresponding part of the MEMF reads, 



£. 



id) 



^9 ^inm-'--' 



mm-'--'m.n 



= {ii'^e^B^, 

+ i"^e,„menei {2eBi^rn - (e + 1) S„,,z) + eff 'enSi,„} (^u<^°^^) 



(28) 
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Figure 2: An illustration of x J effect in disk geometry. Direction of rotation is marked by 
rj, the large-scale toroidal field has opposite direction to rotational velocity and it is marked by 
0, what means that LSMF is perpendicular to the figure's plane and it is directed to the reader. 
The loop of fluctuating magnetic field, b', comprises LSMF that is nonuniform along the axis of 
rotation. Its direction is marked by double arrows. The small-scale Lorentz forces induce the 
azimuthal fluctuations of velocity, u' ~ (b'-V)B. They arc marked by dashed lines ending with 
arrows. The Coriolis force deflects these fluctuations to radial direction (this is marked by dotted 
lines). The resulting electromotive force has the same direction as the original LSMF and it is 
proportional to (6'^) (f] • V)B. 



where functions / 



(d) 
{"} 



(fl*) are given in Appendix A. If we put the magnetic fluctua- 



tio ns in background turbu lence equal to zero in ([28]) (e — 0), we return to results obtained 
by Kichatmov et al.l (1994). The magnetic fiuctuati o ii contributions in (pS]) give rise to the S7 x J 
effect fsee lRadleil (|l969f ): lKrause and Radleil (|f980l ): iRadler et ahl (|2003l ): lKichatinovl (|2003l )) and 
to additions in anisotropic diffusion. In the slow-rotation limit eq. ([28]) can be reduced to 



e„ ((e + 5) Bna + QeBi^n] 



n* 
To 



(29) 



Eq. ((29]) corresponds to results bv iBrandenburg and Subramanian (20o3)- Note, only magnetic 
fluctuations contribute to the induction term (e-V) B. The physical interpretation of this effect is 
straightforward and it is shown on Fig|2] 

Lets consider the situation in disk geometry and the rotating media penetrated by the inho- 
mogeneous toroidal LSMF. For simplicity, we assume that LSMF is nonuniform along the axis of 
rotation. Let the direction of LSMF be opposite to direction of rotating plasma. If the loop of 
the small-scale fluctuating magnetic field comprises LSMF, it induces fluctuation of velocity in az- 
imuthal direction. The influence of the Coriolis force declines the velocities in radial direction. The 
effective electromotive force is co-lined with original LSMF and is proportional to (5^) (SI • V) B, 
see FigH 

The shear-current effect Rogachevskii and Kleeorinl ( 20031 ) (hereafter RK03) is of similar nature 
because the large-scale vorticity W = V x V and the Coriolis force act on the turbulent motions 
in a like manner. The additional contributions due to shear in the diffusion part of the mean 
electromotive force are expressed as follows. 



7{V) 



,(0)2 



(30) 



where Ci = (e - 3/5) r^/G, C2 = (e - I) t^/?>,C3 = (! + £) r^/lS,^ = -(?£ + II) r^/SO. Coeffi- 
cients Ci_4 correspond to those from RK06. After applying the mixing-length approximation to 
RKOe's results we get Ci = -2t^/5, C2 = -4t^/15,C3 = 0,6*4 = -t^/^. In confronting these 
coefficients to ours, we see the difference. It can be explained, in part, by the crudeness of the given 
variant of tau approximati o n. Th e comparison with RS06 is given in Appendix B. As shown by 



Rogachevskii and Kleeorinl (|2007l ) the spectral t approximation is capable to give result in closer 
agreement with those of SOCA. 
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In the commonly accepted scheme of the solar afl dynamo, the poloidal LSMF of the Sun 
is produced from the large-scale toroidal magnetic field via the alpha effect. Expressions (|28|30|) 
hold contributions which are capable to induce the MEMF along the LSMF and consequently 
these terms are potentially very important for the solar dynamo because they provides additional 
induction sources of the large-scale poloidal magnetic field of the Sun. Below, I consider the 
efficiency of induction effect along the nonuniform LSMF due to global rotation and shear. 

For the sake of simplicity we restrict consideration to the axisymmetric LSMF in the Keplerian 
disk in the disk geometry. In cylindrical coordinates (r, 0, z) the axisymmetric LSMF can be 
expressed via B = iJe^ -I- rot (^e^) and the global rotation velocity is V = rfte^. We assume 
that toroidal LSMF exceeds its poloidal counterpart, B « iJe^. In (|28l30p we leave only those 
terms that induce the toroidal MEMF and skip the usual contributions due to turbulent diffusion 
as well. 

In the Keplerian disk we have 9 log $1/9 log r = —3/2. For the given conditions the contribution 
of shear in [3D] is defined by terms at 6*3,(74. It is calculated as follows 

« (C3(V,V)4V,B)^ + C4(V,B)^(V^F),) , (31) 

where covariant derivatives are (VrV)^ — rdr{V ^/r),{y ^Y)r ~ —V^jr and (V^B)^ = dz^- 
Then, the contribution of shear to the MEMF is defined by rdr {V ^^rji'') Tc = — 3il*/4 and 
—TcVtfj/r — —.hVl*. Our derivations are valid in the case of the weak shear flow, |Vi,jrc| <C 1. 
For the Keplerian disks this condition is fulfilled if 17* <C 1. In taking the latter into account 
and using we find the azimuthal component of the MEMF generated from the non-uniform 
toroidal component of LSMF via effects of the global rotation and shear, 

5,«^(2£-M)(.(°)^)r.g. (32) 

Therefore if the LSMF is concentrated to the plan of disk the induced MEMF is in direction of 
the LSMF. 

3.2 Slow rotation, arbitrary LSMF 
3.2.1 Spatially uniform LSMF 

In this part of the paper we consider results obtained for the slow rotation limit. In what follows, 
no restriction is applied to the strength of the LSMF. The MEMF, that is induced due to influence 
of rotation and stratification on the turbulence, is described with expression 

= (^.(°)2)r,{^('^)(GxB)+^('^)(UxB)+r,(f|.B) (^i")G + ^("o^u) (33) 
+ TeB ((^^."^ {n-G) + (fl-U)) + (^^^^^ (B • G) -I- (B • U) 



""^^ (^t^ (B . G) + (B . U)) + ..4^('^)R 



B^ 



wh ere are functions of 8 defin e d in the appendix. This f o rmula generalizes the similar results 



bv iRiidiger and Kichatinov (1993); Kichatinov and Riidiger ( 1992 ) taking the density stratifica- 



tion, magnetic fluctuations and cur rent helicity into account. The no nlinear MEMF of helical 
MHD turbulence was considered bv iRogachevskii and Kleeorin ( 2004a ). as well. For the strong 
LSMF limit we obtain 

^ {i((. + i)(B.u) + i(5^(B.c))(a-(5|)5) (34, 
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cp^"' (£=1) 
RK1993 



Figure 3: The quenching functions for isotropic components of a-effect. 



The results by iRiidiger and Kichatinov (Il993h can be recovered from (j34p . if we put G = and 
e = 0. Following to arguments given in the paper cited above, we conclude that the MEMF like 
does not produce a dynamo. 



Th e first term at the upper line of (f33|) describes the so-called "turbulent buoyancy" (jKichatinov and Riidigei 
The expression (|34p shows that the transport of LSMF is downward for the strong mag- 



netic field limit. For the case of the weak field we get t/?^"'' ~ e/6 + (6e — 8/3^) /15. Then, if we 
neglect contributions due to small-scale magnetic fluctuations, we obtain that for the weak field 
transport is upward (opposite t o direction of G) . In this case the effective drift velocity is propor- 
tional to the LSMF's pressure ( Kichatinov and Riidige ^ (Il992h ). In this aspect it is similar to the 
usual buoyancy of magnetic flux tubes (jParkerl (|l979f l~ Furthermore, we find that the large-scale 
inhomogeneity of magnetic fluctuations provide the downward drift of LSMF in the whole range 
of magnetic field strength. 

The quenching functions for the isotropic components of a effect are shown on Fig|31 There, for 
comparison, via the dash-dotted line, we show th e curve corresponding to quench ing of isotropic 
components of a effect obtained within SOCA in ( Riidiger and Kichatinov ( 19931 )). 

In the strong LSMF limit we found th at a effect is quenched as which is different from 
results by iRtidiger and Kichatinov ( 19931 ) and similar to findings by iRogachevskii and Kleeorin 
(l2004al) . Though, as seen from the figure, the numerical difference between the quenching curves 
obtained within SOCA (dash-dotted line) and MTA (dashed line) is within a few percents. 
The non-linear electromotive force induced by shear is expressed as follows. 



^1 



is) BiBk-. 



B 



BlBr 



■{Vtn~Vn.l) (35) 



+ T^/i^") \ + ^f'^ (F„,, + F,,„)S„ + ^^"^ (W X B) 

B 



From the structure of (j35p we can conclude that contributions with c/?2*3 5 and the second term 

in brackets with ip^^^ can be interpreted as the a effect. The terms with \ and provide 
the pumping of LSMF. Surprisingly, the a-effect like terms survive even in the limit of the strong 



11 



0.15 
0.1 
0.05 


-0.05 
-0.1 

-0.15 
-0.2 

-0.25 
-0.3 

-0.35 








"" ■ 




/ ^ ' 


, - ' 

93 — 






CPi (£=0) 






(Pi(e=l) - 






92(6=0) ■ 









0.1 



10 



Figure 4: Functions defining the nonlinear turbulent diffusion of LSMF (see eq ([33 
magnetic field. In this case we get 

^'t^\si^'x = £inm ^^^^ ^ !_2^ Vl,kUnBm — BiVl^rnUn 

+ ie+ 1)(U . B)^(F,„ - F„,0} (u^°'')r? 



(36) 



2, (0) /-FT B„iB n 



T.-/ir <! ^^™.n^J^B^ - (V^n.^ + V,,n)Bn + (W X B)^ }> , 



According to ([27]) and ([36]) the pumping of the LSMF due to joint effect of current helicity and 
shear have the same sign for the weak and strong LSMF. 

3.2.2 Diffusion, fi x J and shear current effect 

The re sults for nonlinear turbulent diffusion arc similar to those found within SOCA bv lKichatinov et "al 
(|l994h . We have 



= |^3V X B + xB^)xB) ^ ^ ^ _| ^ ^(^^^ ^^^^ 

where f stands for the contributions due to rotation. The corresponding quenching functions 
are given on FiglH 

The next formula generalizes the results for the nonlinear diffusion of LSMF to the case of the 
slowly rotating media, 



^^'"V, (f2 • B) + v^i^ 



{w) 



(n -B) 



+ ^^^^ 



B-V 



log 



(B> 



2 

(to)- 



-vaog (s'j + 's, 

(r2 • B) (B • V 



(B • V) • B) 



■^3 ^» -2 2 



<^6 ' -^^JT^ (B ■ V) + <^^'"^B,^^^ log 
B ^ 



log (s') (38) 
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The last two terms at the third hne in ([SS]) are related with the generation of MEMF along the 
direction of LSMF. The corresponding functions ip^^^ and {p^"'' are shown on FigO As can be 
seen there, in the absence of the background magnetic fluctuations (e = 0) the generation due to 
VL X J-efl'ect exists only in nonlinear regime. 

If /3 > 1, functions ^^2"^ and (py™"* have opposite signs everywhere. Note, while the term 
(O • V) Bi induces MEMF in direction of LSMF's gradients along axis of rotation, the term 

Bi [Vl ■ V) log {b^^ induces MEMF in opposite direction. Formally, the latter effect is similar 
to a-effect. The only difference with the standard a is that instead stratification parameters of 
turbulence we have a parameter which is related with nonuniform distribution of the LSMF's 
energy. For the solar magnetic fields the effect is antisymmetric about equator. Below, it is shown 
that in the strong LSMF this a is quenched by factor j3-^ which is lesser than for standard a. 
For the limit of the strong LSMF we get 

fi™^|;3-.oo = |(17e + 47)(v.(f^.B)-i^^Vaog(#)) (39) 



/ (B • V) /_2\ _ (B • V) (1] • B) in-B) ,- , _ ^ 
(2l£ + 43) I ^^^^ log [b j + B,^ >- + (B • V) B, 

fl7-B)(B-V) /-2\- 
3 21e + 43 ^ '-A^ log [B ] B, 

2B ^ > 

(37e + 27) fe^^^ log (#) - • V) B 



\ 2 ° V / ' 7 J 512/3 

From there we find that 17 x J-effect maintain the generation part of the MEMF even for the 
strong LSMF. The amplitude of effect tends to constant as the strength of LSMF is increased. It 
hardly possible to make a definite conclusion about the dynamo effect in this case because the 
generation part of (I39p is contributed by terms with opposite signs. 
The MEMF's contributions due to shear are defined by. 
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BlBr, 

f_ 

I ( (^)T7 1 (^)T7 ^"D I (v)BiBnjj -= {v)BkBi— — 

(y-^BkBi— — (y)— — BkBi— (y)— (y) — 

+ 'y^'s _2 V^m,«-Pi,fc + ^9 V i,n-0;,m + _^ (^lo -^fc^m + ttOn -Pm,fc, 

where, for the sake of simpUcity, we leave only the largest contributions and those which are 
important for the solar-type dynamo models, where the strength of LSMF component along direc- 
tion of the large-scale flow dominates components directed along the shear. Reader can find the 
expressions for ip'n'^ in Appendix A. The full expression has a much more complicated tensorial 
structure than (|in)l . In the case of the strong LSMF we get, 

^ri^-^oo = ^f{u^'^^)s.n.^\J^-^Vr.,-'f^l,^|B„,, + {s (41) 

— ;:;7; Vn,lBi,„i ^ur-Vl kBk^m H ttt-Vi nBl m 

zL) z iU 

BkBi— ( , -.N^ ^ + 3^ 

+ _2 ^ l,n \ + ^)^m,k — ^k,m 

Now, we would like to consider efSciency of induction effect along the nonuniform LSMF due 
to global rotation and shear in nonlinear regimes for the Keplerian discs. As before, we assume 
a disc penetrated by the large-scale toroidal magnetic field that is nonuniform along the axis of 
rotation. From (|39l40p we get 

where the quenching function is (^(""^^ = (^2™'' + 'p't"'' + ^-^{Vt^^ + Vw^) + 'y^i^''- Note, eq. (j42|) 
transforms to eg. ([5^ in limit /3 ^ 0. The dependence of (/j^™^) on the LSMF's strength is shown 
on the FigEl 

Results given on the Fig[S] show that the ip^"^^^ is positive for /3 < 1 and negative for /? > 1 
for all e. This result supports an idea about the change of dynamo type in passing from linear 
to non-linear regime of the LSMF's generation by x J and shear-current effects. Previously we 
found that the induction term due to x J effect tends to constant when (3 ^ oo (eq. p9|l ) while 
the induction term due to shear-current effect is growing under (3 oo (eq. (|4ip ). Therefore, 
the primary nonlinear generation effect in the differentially rotating uniform MHD turbulence 
penetrated by the nonuniform toroidal LSMF may be due to shear-cu rrent effect. The same non- 
linear d ependence of shear-current effect was discovered in the paper bv lRogachevskii and Kleeorin 



(|2004ah for the different kind of MTA. In the next section I show that the given sources of the 



MEMF ultimately result to current helicity generation. Therefore, the effect considered above is 
saturated dynamically due to magnetic helicity conservation law. 



4 The current helicity evolution 

As we have seen, the current helici ty contributes to the different kind of M EMF's action, not only 
to the a effect. The recent papers ( Subramanian and Brandenburgl [20041 ) show that the magnetic 



helicity conservation law can be described in terms of the current helicity evolution if the assump- 
tion of the scale separation is fulfilled. For the time being the redistribution of current helicity over 
the space scales is not satisfactory understood. One attempt to describe the helicity evolution in 
turbulent media penetrated by LSMF was given in the papers by iBrandenburg and Subramanian 
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Figure 6: The dependence of induction effect along the nonuniforni toroidal LSMF on the strength 
of magnetic field. 



( 2005 ): Subramanian and Brandenburg ( 2004 ). Here, we will follow their results and obtain the 
explicit evolutionary equation for the current helicity. The equation in question can be derived 
from pi4p . After integration over the large-scale variables we can get the general equation for the 
current helicity in the following form, 



dhc 
dt 


Th 


2 

MP 




+ ifcpVn 


Kin 

V P 



fc^i.,p^ - i>ci, (k . V) (^^) - ^ (k • V) (^>ci,^ ] (43) 



1, 



dk. 



The third order moments were replaced by ~hc /t^ , t^- is a relaxation time for the current helicity 
. This is a rather rough way because the tripl e correlation s may gi ve important contribution for 
the helicity redistribution over the space scales ( Frisch et al. . 1975; Kleeorin and Ruzmaikinl . ll982l : 
iKleeorin and Rogachevskii . 1999f) . Because of the very rough assumptions used in derivation of 

it should be considered with caution. In spite of the latter, the equation provides a useful 

tool f or investigation the nonlinear saturation in helical mean-field dvnamo(Brand enburg and Subramanian 
l2004h . Except for contributions due to density stratification and shear, equation (j43|) can be re- 



produced from results of BS05 after substitution identity eijk£ipqSqim = £imkSjp — simj^kp in eq. 

7^ and use we can write the 



(10.71) there. Inspection of shows that if we replace 
evolutionary equation in the following form. 



dhc 
dt 



2 (f • B) 



Th 



Brn 
P 



PP 



-ii</pfc"V 
1„ 



p 



(44) 



dk. 



According to (jFrisch et ah ^ 1975t Kleeorin and Ruzmaikin , 1982 ; Vainshteinl . Il983 ; Brandenburg , 
2001; Vishniac and Cha 2001 ) the first term in (|44p is responsible for helicity g eneration in turbu- 
lent medium. The rest part of equation can be inter preted as the helicity fluxes (jVishniac and Choi . 



20011 ISubramanian and Brandenburel . I2004L l2005h . The given expression for helicity fiuxes is in 



complete because the contribution of the third order moments is dropped in . As the first step 
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we consider the case of the weak LSMF. From (jH]) and p8|17|19p we get 
^ {fi'^B' + f'f (e • B)^) + 2/f ) (B . (V X B)) 



dhc , 



(45) 




(e.B)(^^/f (B.G 
\f^^ (e . B) (e . (V X B)) - ^^^^-^ (e • B)^ - f^^^^ 



e-B 



where substitution (u^'^'^)^^^ t^"^ was used, and e = r2/|ri| . Here, we dropped the contribu- 
tions due to shear because their effect to the mean electromotive force was computed only to the 
zero order terms about angular velocity. Furthermore, in ()45|) we kept only those contributions 
which could be the most interesting from the stellar dynamo applications standpoint. Note, for 
the equipartition case, e = 1, helicity evolution satisfies the simple equation: 



dhc 
dt 



1 

Th 



-he = - 



(46) 



It is in accordance with equation for the magnetic helicity density obtained bv lSubramanian and Brandenburg 
(fioOS) . As an example of application of (j46| to the problem of the nonlinear saturation of alpha- 
effect, consider the dynamo in the fast rotation limit. For the sake of simplicity we restrict 
ourselves only with the isotropic components of a-effect and neglect the helicity loss due to hc/Th- 
From gSl) and (HH) we get 



dhc 
dt 



4t, 



,(0)2 



Vt* (e • G) - he 



(47) 



where we keep the contributions of order Q,*^^ for the current helicity, and drop the terms which 
are due to nonuniform LSMF. If L is the typical spatial scale of the LSMF then the eq. (|47p 
is justified when LGVL* > 1 and ^J.p\hc\ » |B • (V x B)|. The point to note that in ([T7|) we 



widely used in the literature (iKleeorin and 


■luzmaikinL 19821 Vainshtein anc 


Kitchatinov. 1983: 


Vishniac and Choi 2001; Kleeorin et al.L 2003 


; BrandenburK and Subramanian 


. 20041). With initial 



condition, t = 0, /ic = 0, we write, similar to 



yaiiishtcin (1983), the solution of eq (|T7)) as follows. 



he = 



,(0)2 



G) 1 - exp 



TT 

'4^ 



ITdt 



(48) 



The given solution shows that under t ^ oo we get he 2W (^u^*'^^) (e • G) t^- On this basis, 
and in taking into account (|26p . we can conclude that a-effect will saturates exponentially under 
the increase of the LSMF strength. Furthermore, this conclusion w as con firmed with numerical 
dynamo model which is considered by author in the separate paper (|Pipinl . l2007 ). 

Next, we consider the equation for the current helicity evolution for the slow rotation limit. 
No restriction is applied to the strength of LSMF. The contribution of shear to the transport 
and generation part of equation is described with a quite bulky tensor expressions and we decide 
to restrict ourselves with terms which have either a finite limit under /3 — > or the amplitude 
functions that are greater than 0.1. We write the evolutionary equation for the current helicity as 
follows: 



dh 



. -he = -^(f •B)+V.i^^Fp,™/ic + (^2G + ^3U).w(^(")2 
dt Th upli B ^ 

+ • ([V'sV X B + V4 (U X B)] (B • V) + VeWS') + (e - 1) {...} , 



(49) 
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Figure 7: The isolines of the angular velocity distribution (left) and the corresponding vector field 
of the large-scale vorticity (right). 

where W — V x V . Quenching functions tp{n} ^re given in Appendix A. Symbol {...} denotes 
those terms which are not important in the case e = 1. Taking the Taylor expansion of for 

the case B ^ (keeping terms) we get 




'p,mhc- ' g ' {^u^^>'j~V ■J' (50) 

([VxB-(UxB)] (B-V)), (51) 

where we apply the equipartition condition, e = 1, as well. The direction of the helicity flux 
due to the first contribution in ([5T|) . = (J^u''^^^'^ /6 + 2i?^/ (15/ip)^ W, depends on distri- 
bution of the large-scale vorticity solely. The second term depends on details of the dynamo 
action. To estimate the direction of the helicity transport due to on the Sun we com- 
pute the vector field of the large-scale vorticity W. In the spherical coordinate system we have 
W = sm0e^dfl/d9 — rs'm9e^dfl/dr, where r, 6 are the radial distance and the pol ar angle, respec- 



tively. The distribution of the angular velocity is taken as an analytical fit given bv lBelvedere et al 



(2000). It is shown at the left side Fig[7|The computed vector field of the large-scale vorticity is 



shown at the right side FiglT] 



The given figure shows the possibility of the outward helicity flux from the dynamo region due 
to shear. Note that one component of the helicity flux is due to the small-scale dynamo, 
(u'-°^^'^'W/6, and another is due to the LSMF, 2B^W / (15^p). Among two the contribution of 

the small-scale dynamo is likely to be dominated in the depth of convection zone. While the flux 
due to the LSMF may be important at near the surface level. T he latter effect may produce the 
significant outward flux of helicity only with the open boundaries (jBrandenburg and Subramanian , 
[2004; Subr amanian and Brandenburgl 2005f ). At the near surface level the amplitude of the large- 
scale vorticity, |W| w 4 x 10~^s^^ sa 1.5 x 10~^da y~^ . The magnitude of the surface magnetic 
flux change during the solar cycle is about 10^'^Ma; ( Schrijver and Harvevl . 1983) ■ Then the mag- 
nitude of the helicity outflow f rom 2 B W/ (15^/9) is about 10'^^ Mx'^day~^ . It is compatible with 
estimations given by iDeVord (|2000h . We have estimated only one part of the helicity flux. The 
numerical dynamo model based on the given results would help to get a more definite conclusions 
about this subject. 
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5 Summary 



In the paper the mean electromotive force of turbulent flows and magnetic fields is computed 
analytically using the framework of mean-field dynamo theory and MTA (minimal r approxima- 
tion). There is an overlap in results obtained with SOCA and with MTA approximation. The 
two approximations give qualitatively the same results about nonlinear dependence of mean elec- 
tromotive force on the strength of LSMF or on the Coriolis number. Also there is a difference 
between predictions of SOCA and MTA for mean-electromotive force expressions if the shear is 
taken into account. This difference can be explained, in part, by the crudeness of the given version 
of tau approximation. The deficient accuracy of calculations of shear contribution is due to an 
assumption about the scale-independent r. In whole, the accuracy of the theory presented in the 
paper is comparable with the mixing-length approximation. It has no the firm grounds and should 
be considered with caution. 

Finally, I would like to focus on the new findings of the paper. In this study it is shown that 
the new interesting component of transport of LSMF appears due to joint contribution of current 
helicity and shear. The effect does not disappear in the strong LSMF limit, (3^1. It may be 
important near the base of solar CZ where the influence of rotation and shear on the turbulence 
is quite strong. Furthermore, the analysis, which we carried out for the current helicity evolution, 
suggests that the shear and rotation may redistribute the helicity in solar CZ amplifying it (in 
amplitude) at the near equatorial regions in agreement with observations. Beside, the effect of 
rotation and stratiflcation on the he evolution is calculated explicitly. Basically, the equation for 
current helicity is obtained using the same approach as for the mean electromotive force and on 
the base of quantities which are explicitly gauge invariants. Therefore, we can expect that the 
dynamo model based on the above approach could be capable for meeting the requirements of 
both solar and stellar dynamo simulations. 
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Appendix A. 

This part of appendix gives the functions of the Coriohs number defining the dependence of the 
turbulent transport generation and diffusivities on the angular velocity. 
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The dependence of turbulent diffusivities on the Coriolis number (eq. ()28p ) is given by 
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The magnetic quenching functions of the generation and transport effects in eq. are 
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The nonlinear turbulent diffusion of the LSMF in ([571 ) is expressed with help of the following 
functions 
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The effect of slow rotation and nonuniform LSMF on the MEMF (eq. is expressed with help 
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The magnetic quenching functions for the shear-current effects are 
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The magnetic quenching functions for the current helicity evolution equation: 
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Appendix B. Comparison with some of results given in the paper by 
Radler & Stepanov (2006). 

This part of the article contains the comparison some of our results with those from RS06. We ap- 
ply the mixing- length (MLT) approximation to expressions obtained in RS06. In this procedure we 
replace the spectrum of turbulent fields by the single-scaled function of the for m 5 (k — S (uj) , 



where £c is the correlation length of the turbulence and put r/k^ = vk^ = ^. (jKichatinovl . Il99l[ ) 



The effect of stratification and shear. The structure of the electromotive force obtained by 
RS06 can be reproduced if we decompose the gradient of the large-scale flow Vij into symmetric 
and antisymmetric parts via 

1 

2' 



V^j = - T^'^ijnWn, (52) 



where Wi — einmVm,n is the large-scale vorticity and Dij — {Vij + Vj^i) /2 is the rate of strain 
tensor. After substitution (152p to ([?7|) we obtain 

^i'' = ( e„,mUkBnD,nkA4 - 4^ (W • B) C/, + (U • B) f A3 - ^ + 4^ + ^) (53) 
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Using ((551) we find 

r-^aP = \r-\A,~A,)^0, (54) 

r.-^^r = -\r-'[A, + \{A,-A,)^=-^^, (55) 

r-^^) - -ir-^ (^A3 + A4 + i(Ai-A2)) =0, (56) 

rc^^""^ - ir-2(3A4-Ai-A2) = -^, (57) 

= ir-2(A4-Ai~A2)^-^, (58) 

(59) 

where we put e = 0. After applying the MLT to results obtained in RS06 we find r~^a$^^'' = 
19/120, = -7/240, r-^^lW') = -1/48, t-^^^'') = -39/120, r.-^gC^) = -21/120. 

The efTect of nonuniform LSMF and shear. For the shear-current effect, after substitution 
of ((52|) to (j30p we arrive to the following representation of e\^'' , 

) = (w . V) + ^4^V, (W.B)} ( (60) 

+ e,:,™ { (Ci + C2) S„,i + (C3 + C4) A„/ 
Using this formula we obtain 



.-2 



= ir-2(C3-C4-Ci+C2) = ^, (61) 

r-'^f'^^ - ^r-2(C4 + C2-Ci-C3)--^, (62) 

T-''^^''^ = -^r-2(Ci+C2 + C3 + C4) = ^, (63) 

T-'/3(^) = -ir-2(Ci+C2-C3-C4)-0, (64) 

where we put e = in Ci_4. After applying the MLT to results in RS06 we find = 1/12, 
^-2~(W) ^ „i/3o^ r-2K(^) = 13/30 and r^^/jC^) = 7/6O. 
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